Structural and electronic properties of Pbi.^Cd^Te and Pbi_ z Mn z Te ternary alloys 
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A systematic theoretical study of two PbTe-based ternary alloys, Pbi-^Cd^Te and Pbi-ajMn^Te, 
is reported. First, using ab initio methods we study the stability of the crystal structure of CdTe - 
PbTe solid solutions, to predict the composition for which rock-salt structure of PbTe changes into 
zinc-blende structure of CdTe. The dependence of the lattice parameter on Cd (Mn) content x in the 
mixed crystals is studied by the same methods. The obtained decrease of the lattice constant with x 
agrees with what is observed in both alloys. The band structures of PbTe-based ternary compounds 
are calculated within a tight-binding approach. To describe correctly the constituent materials 
new tight-binding parameterizations for PbTe and MnTe bulk crystals as well as a tight-binding 
description of rock-salt CdTe are proposed. For both studied ternary alloys, the calculated band 
gap in the L point increases with x, in qualitative agreement with photoluminescence measurements 
in the infrared. The results show also that in p-type Pbi-xCd^Te and Pbi-^MnrTe mixed crystals 
an enhancement of thermoelectrical power can be expected. 

PACS numbers: 71., 71.15.-m, 71.20.Nr, 71.28.+d 
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I. INTRODUCTION 

Since 1959, when the first Radioisotope Thermoelec- 
tric Generator (RTG) was presented, PbTe attracts 
constantly a lot of interest due to its thermoelectric 
properties f^t— PbTe is also widely used for mid-infrared 
lasers and detectors. In PbTe the narrow direct band 
gap at the L point increases from 0.19 eV at 4.2 K to 
the value of E g = 0.31 cV at room temperatures^ This 
allows excellent band structure engineering. Recently, 
also ternary systems based on lead telluride are of con- 
siderable scientific interest because of their potential in 
device applications £ An important characteristic of the 
PbTe-based ternary alloys is that their band gap is very 
sensitive not only to temperature, like in PbTe, but also 
to composition. It has been shown that the energy gap 
of PbTe increases monotonically when alloyed with Cd&£ 
as well as with Mn£ 

Due to big, ca 1 eV, difference between the energy gaps 
of PbTe and CdTe, PbTe-CdTe system appears especially 
suitable for band-gap engineering with a potential for 
a variety of photonic, thermoelectric and photovoltaic 
applications^ The growth of uniform Yhi_ x Cc\ x Te sin- 
gle crystals is, however, very limited by extremely low 
mutual solubility of both materials^ The latter results 
from the difference in the crystal structures - lead tellu- 
ride crystallizes in rock-salt structure while cadmium tel- 
luride in zinc-blende structure. The limited mutual solu- 
bility of PbTe and CdTe was yet exploited for obtaining 
PbTe quantum dots in a CdTe matrix. This was achieved 
by thermal annealing of two-dimensional PbTe epilayers 
embedded in CdTcAi In Ref. [Tl] it was shown that the 
size of the dots can be controlled, what allows for tuning 
of the quantum dot luminescence over a wide spectral 
range. As a result, ultrabroadband emission from a mul- 
tilayered quantum dot stack was demonstrated, which is 
a precondition for the dcvclopnrcnt of supcrluminescent 
diodes operating in the near infrared and midinfrared. 



On the other hand, it is expected that thermoelectric 
properties of PbTe should be improved by implement- 
ing CdTe nanostructures in the material. In Ref. [l2| a 
possibility to increase the thermoelectric figure of merit 
parameter ZT of certain materials by preparing them in 
the form of quantum- well superlattice structures was pre- 
dicted. It has been also shown that CdTe nano-clusters 
embedded in PbTe lead to considerable changes of the 
derivative of the carrier density of states at the Fermi 
level and can influence the thermoelectrical properties of 
the material^ These theoretical results together with 
the recently reported fabrication of CdTe quantum dots 
in a PbTe matrix^ open doors for using PbTe-CdTe 
structures to enhance the performance of thermoelectric 
devices. 

Despite the mentioned above difficulties, bulk 
Pbi_ 2 Cda;Te solid solutions in the form of polycrys- 
tallinc samples were obtained by both, the Bridg- 
man technique^ and by a rapid quenching fol- 
lowed by annealing ,SiZii£ Recently, high quality sin- 
gle Pbi-aCd^Te crystals with x as high as 0.11 were 
obtained^ by self-selecting vapor growth method^ 
These efforts were motivated by one more advantage in 
using PbTe as a base for forming ternary alloys, i.e., by 
the fact that in these materials the relative contributions 
of light and heavy holes, thus the electrical and opti- 
cal properties of the system, can be tuned by changing 
the composition or temperature. Indeed, it was shown 
that while the energy gap of PbTe increases with Cd 
content, the energy separation between the light and 
heavy hole valence bands is considerable rcducedi&i Sim- 
ilar behavior was observed also in Pbi-^Mn^Te crystals. 
In contrast to PbTe-CdTe system, the solid solution of 
PbTe and MnTe leads relatively easy to Pbi-^Mn^Te 
single crystals with x up to 0.10. The experimental 
studies of PbMnTe suggest that in this material adding 
Mn ions to PbTe changes also the relative positions of 
different valence band maxima, offering a possibility of 



improving the thermoelectric properties. In the p-type 
Pb 1 _ :E Mn : j.Te crystals it was shown that at room tem- 
perature and for constant carrier concentration (p — 
2 x 10 18 cm -3 ) the thermoelectric power increases rapidly 
with the increase of Mn content, thus improving the ther- 
moelectric figure of merit parameter Z.— 

In this paper we present a systematic study of the 
structural and electronic properties of Pbi-^Cd^Te and 
Pb 1 _ ;c Mn : j.Te ternary alloys, for which either ab initio 
or tight-binding methods were used, when appropriate. 
The first principle calculations of the stability of mixed 
crystals are presented in Sec. II. In Sec. Ill, improved 
tight-binding description of the valence and conduction 
bands of PbTc, rock-salt CdTe and MnTe bulk crystals 
and the results obtained for the band structure of their 
solid solutions are shown. Section IV contains our conclu- 
sions, in particular the predicted, at using the obtained 
band structures, thermoelectric properties of the mixed 
crystals are discussed. 



II. STABILITY OF THE MIXED-CRYSTALS 

In Pbi-^CdzTe mixed crystals the transition from a 
ten-electron (x—0) to an eight-electron system (x=l) oc- 
curs, which is accompanied by a change in the crystal 
structure from rock-salt (RS) to zinc-blende (ZB). The 
fundamental question to be asked for these crystals is 
the x value of this structure change. To answer this 
question, we analyze the stability of the crystal struc- 
ture in Pbi-xCd^Tc alloys using ab initio density func- 
tional theory (DFT) method. Moreover, using the same 
method, we determine the dependence of the lattice pa- 
rameter on Cd and Mn content x in Pbi_ x Cd x Te and 
Pbi-^Mn^Te crystals, respectively. The calculations are 
performed within the Vienna ab initio simulation pack- 
age (VASP)i 19 ' 20 For the atomic cores the projector aug- 
mented wave (PAW) pseudopotentials^i are used. The 
exchange correlation energy is calculated using the local 
density approximation (LDA). The atomic coordinates 
are relaxed with a conjugate gradient technique. The cri- 
terion that the maximum force is smaller than 0.01 eV/A 
is used to determine equilibrium configurations. In all 
the calculations, the energy cutoff is set to 16 Ry for the 
plane- wave basis, which is sufficient to obtain converged 
structural properties. Since the impact of nonscalar rela- 
tivistic effects on the structural features is negligibl o 22 i 23 
we do not take these effects into account. Properties of 
Pbi_ x Cda,Te and Pbi-^Mn^Te systems are calculated 
with (2x2x2) simple cubic supercells containing 64 
atoms. The Brillouin zone integrations are performed 
using 4x4x4 Monkhorst-Pack k-points meshes. All 
the atomic positions and the volume of the supercells are 
calculated with relaxation and re-bonding allowed. Our 
study is carried out for zero pressure and zero tempera- 
ture. 

To study the stability of the crystal structure of 
Pbi_ x Cda,Te alloys we calculate the total energy of RS 



PbTe and ZB CdTe supercells, in which we exchange 
successively the Pb(Cd) cations by Cd(Pb) ions, respec- 
tively. Thus, in the supercell of RS structure for Cd 
content x = we have a PbTe crystal while for x = 1 
we obtain a RS CdTe crystal. On the other hand, con- 
sidering the ZB supercell, for y = (1 — x) = we obtain 
the CdTe crystal and for y = 1 a hypothetical ZB PbTe 
crystal. One Cd atom in the RS PbTe supercell cor- 
responds approximately to the Cd content x = 0.03 in 
RS Pbi-^Cd^Te alloy (similarly one Pb atom in the ZB 
CdTe supercell leads to ZB Pbo.03Cdo.97Te). 



-5.5 



> 

— -6.0 



TO 



-6.5 



-7.0 



-7.5 



8 

E 
c 

~^ 
o 

llT 



-8.5 



-9.0 



- 


ie ^X 


* -kr^ 


* RSPb^CdTe 
» ZBPbCd^Te 

11 ' bulk bulk 







0.0 



0.2 



0.4 0.6 

Cd content (x) 



0.8 



1.0 



FIG. 1. (Color online) Dependence of the total energy of 
RS Pbi-^Cd^Te (black stars) and the ZB Pbj,Cdi_„Te (red 
dots) on the Cd content. The total energy is normalized to 
the number of cation-anion pairs. This energy is compared 
with an appropriate for the Cd content x sum of the energies 
of two separate, RS PbTe and ZB CdTe, phases (blue line). 

Fig. [H presents the dependence of the total energy 
(taken per atomic pair) on the Cd content for RS and 
ZB Pbi_ x Cd x Te crystals. First, we see that the energy 
difference between RS and ZB structure is much bigger 
for x — than for x = 1 (the energy of ZB PbTc is higher 
than that of RS PbTe by ~ 0.666 eV/atomic pair, while 
for CdTe the RS structure leads to higher energy than 
ZB by ~ 0.194 eV/atomic pair). We also observe that 
the RS structure of PbTe is preserved in the calculations 
for Pbi-zCd^Te mixed crystals practically in the whole 
range of Cd concentrations. In contrast, adding Pb to 
CdTe destroys very quickly the ZB structure - after re- 
laxation the structure is maintained only for Pb content 
up to ca y = 1 — x = 0.25 (in our PbCdTe supercell 
with 32 cations, this corresponds to 8 atoms of Pb and 
24 atoms of Cd). Therefore, in Fig. Q] only the points 
(red dots) which correspond to the preserved ZB struc- 
ture of Cdi-yPbyTc crystals are shown. These results 
are consistent with the observation that it is easier to 
obtain Pbi-^Cd^Te crystals than Cdi_ y Pb y Te. Indeed, 
while there are no reports on CdTe highly doped by Pb, 
the successful growth of Pbi-^Cd^Te monocrystals with 
Cd content up to x — 0.11 was reported^ 

We want to check for which Cd concentration the tran- 
sition from RS to ZB structure should occur. As one can 



notice in Fig. [T] the total energy of the mixed crystal 
for RS structure is lower than for ZB for nearly whole 
range of Cd concentrations, i.e., up to i ~ 0.8. This 
means that virtual Pbi^Cd^Te crystals with such high 
Cd concentrations would still have RS structure. How- 
ever, for any x value this energy is slightly higher than 
the total energy of separate phases (blue line in the Fig. 
HJ which we calculate as: (1 - x)E^ e + xE^ e . Still, 
for x up to ~ 0.25 the energy difference between the line 
representing the energy of separated phases and energy 
of the mixed crystal in RS structure is lower than ksT 
in the growth conditions (~ 0.1 eV). This result may de- 
note that the highest Cd content in samples presented in 
Ref. Il3 is close to a fundamental solubility limit. 

In the next step we study the dependence of the lattice 
parameters an of Pbi-^Cd^Te and Pbi-^Mn^Te crys- 
tals on the composition x. For this purpose we con- 
sider a PbTe supercell in its stable RS crystal struc- 
ture, in which we successively replace the Pb atoms by 
either Cd or Mn ions. The dependence of the calculated 
lattice parameter of Pbi-^Cd^Te on x is presented in 
Fig. [5] and for Pbi.^Mn^Te alloys in Fig. [3] We see 
that the relation between ao and x < 0.12 for both 
mixed crystals is linear. The lattice parameters dimin- 
ish with x like da a /dx = -0.434 A (Pbi-sCdxTe) and 
da a /dx = -0.683 A (Pbi_ x Mn x Te). It is well known 
that the DFT calculations underestimate the lattice con- 
stants - here the absolute values of the calculated lattice 
parameters are for both materials about 0.1 A lower than 
the measured. Still, the obtained theoretically rate of de- 
crease of the lattice parameter with x for both crystals is 
in very good agreement with experimental data (compare 
Ref. El for Pbi-zCd^Te and Ref. [H for Pbi-^Mn^Te). 
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FIG. 2. (Color online) Lattice parameter of Pbi-^CdzTe 
mixed crystal as a function of the Cd content x. Green and 
blue dots denote the calculated and experimental values, re- 
spectively. 







Pb, Mn Te 


o 

m 


6.44 


^^-^^a a Theory 

^^-^^^ 9 Experiment - 


it 




da rf 'dx = -0.632 A ^*^^^^ 


■* — ' 
£ 

(0 


6.40 
6.36 


*^\ 


<1> 
O 


6.32 


^^-«^^ 


_l 


6.28 



da^dx = -0.683 A ^^i 




00 0.02 0.04 0.06 0.08 0.10 0.12 






Mn content (x) 



FIG. 3. (Color online) Dependence of the lattice parameter 
of Pbi-^Mn^Te on the Mn content x. Black and red dots 
denote the calculated and experimental values, respectively. 



III. BAND STRUCTURE OF THE 

CONSTITUENT MATERIALS AND THEIR 

ALLOYS 

As shown above, the DFT calculations give much too 
low values of the lattice parameters. It is also well known 
that both, the local density and generalized gradient 
(LDA and GGA) approximations to the DFT method, 
widely used in the band structure calculations, underes- 
timate the fundamental band gaps. In the case of narrow 
gap semiconductors, like PbTe and other lead chalco- 
genides, this leads even to a change of the sign of the 
gap, i.e., to an erroneous "inverted" order of bands at 
the L point of the Brillouin zone. Several attempts to 
overcome this drawback of the ab initio calculations, by 
using GW method 2 ^, hybrid functional^ or by simple 
rigid upward shift of the conduction band^, were re- 
ported in the literature. Recently, it was shown 2 ^ that 
correct band structure of PbTe and correct changes of 
the electronic structure of PbMnTe with the concentra- 
tion of Mn ions and pressure can be obtained within DFT 
by reducing the spin-orbit strength for Pb Qp electrons 
by approximately 40%. 

Here virtual-crystal approximation and tight-binding 
method are exploited in the analysis of the electro- 
nic structure of Pbi-^Cd^Te and Pbi-^Mn^Te alloys. 
Tight-binding methods have proven to be very useful in 
studying the electronic properties of solids. In the em- 
pirical tight-binding Hamiltonian, the matrix elements 
between orbitals centered on different sites are treated 
as parameters, which are adjusted to any known, often 
experimental values. The band structures of the IV- 
VI semiconductor compounds obtained by tight-binding 
metho d 29 ' 30 are widely used to explain the observed phe- 
nomena in these materials. To describe the band struc- 
ture of the PbTe-CdTe and PbTe-MnTe solid solutions 
we need the tight-binding paramctrization for all the 
constituent compounds, i.e., PbTe, CdTe and MnTe. 



We recall that the typical structure of Pbi-^Cd^Te and 
Pb 1 _ rr .Mn :E Te mixed crystals is the RS phase. Therefore, 
to examine the properties of these systems, we need pa- 
rameters describing not only lead telluride, but also cad- 
mium telluride and manganese telluride semiconductors 
in RS structure. While PbTe crystalizes in RS structure 
and RS MnTe is also common, CdTc changes its struc- 
ture from ZB to RS only under high pressure. Therefore, 
the tight-binding paramctrization for the two former ma- 
terials can be found in literature, but to our best knowl- 
edge, there is no tight-binding parameters available for 
RS CdTe. 



A. PbTe 

Let us start from the band structure of PbTe. A care- 
ful analysis of the tight-binding parameters available in 
the literatur o 29 i 30 shows that although they recover cor- 
rectly the band structure, they do not lead to effec- 
tive masses, which are determined experimentally. Thus, 
in order to fit our tight-binding model to all existing 
experimental results, we have performed a new tight- 
binding paramctrization for the PbTe crystal. To de- 
scribe PbTe material we use the sp 3 atomic orbitals, with 
the spin-orbit coupling included. In our model we con- 
sider the nearest-neighbor cation-anion as well as next 
nearest anion-anion and cation-cation intcr-atomic cou- 
plings. We use the experimentally determined energy 
gap in low temperatures, i.e., E g = 0.19 eV^ and as- 
sume that at T = OK the second valence band maxi- 
mum along the S-line is located about 0.17 eV below 
the top of the valence band at the L point, as suggested 
in Ref. l3ll Another experimental input to our fitting 
comes from Ref. |32L where the values of the longitudinal 
and perpendicular PbTe effective masses at the L point 
of the Brillouin zone are given. Spin-orbit parameters 
are matched to the atomic values for Pb and Te atoms, 
which are equal 1.273 eV and 0.840 eV, respectively^ In 
our fitting procedure we allow for an adjustment of the 
latter parameters, but we keep the ratio of cation to an- 
ion spin-orbit coupling parameters equal to the Pb/Tc 
atomic values rate, i.e., to 1.51. Our final values for 
spin-orbit parameters differ from the atomic spin-orbit 
couplings by less than 10% (compare: Table IJ). In Ta- 
ble U the obtained in this work tight-binding parame- 
ters for all, PbTe, RS MnTe and RS CdTe, materials 
are shown in the standard Slatcr-Kostcr notation^ The 
band structure of PbTe resulting from our tight-binding 
paramctrization is presented in Fig. 01 As one can see 
in the Figure, our improved paramctrization gives the 
correct value of the energy gap at the L point and the 
appropriate energy position of the second valence band 
maximum at the S. Moreover, in contrast to the earlier 
calculations, our model provides proper values of the lon- 
gitudinal and perpendicular effective masses. In Table [TT1 
the comparison of the effective masses obtained within 
our model and the other theoretical approaches with the 



TABLE I. Nearest- and next-nearest neighbors tight-binding 
parameters and the on-site energies for RS PbTe, CdTe and 
MnTe crystals. All energies are given in eV and the energy 
zero is always assumed at the top of the valence-band. The 
A denotes the spin-orbit coupling parameter. 



Parameters (in eV) 



PbTe 



CdTe 



E Pc 
Ed c 
E Sa 

Ep a 
ScSaCT 

s c p a a 

PcSaO 
PePaCT 
PcPaK 
dcSaCT 
dcPaO 
dcPalT 

s c s c a 

ScPcO 

p c s c a 

PcPcV 

PcPcK 

dcScCT 
d c p c a 
dcPcK 
d c d c a 
d c d c -K 

s a s a a 

SaPaV 
PaSaV 
PaPaCT 
PaPaK 

A c /3 
A a /3 



MnTe 



8.6528 


-2.3095 


1.6517 


1.2711 


3.7778 


4.5995 




-8.0291 


0.1676 


9.4379 


-9.7514 


-9.5918 


0.8324 


1.0687 


0.2887 


0.3265 


0.9943 


0.9823 


0.0838 


1.3711 


1.9093 


0.2148 


0.9478 


-0.0786 


1.6702 


2.0861 


2.5137 


0.1149 


-0.6885 


-0.2744 




0.2070 


0.1959 




0.7366 


0.5445 




-0.0709 


0.5420 


0.2444 


-0.0951 


0.0954 


0.4909 


0.2331 


-0.1853 


0.4909 


-0.2331 


0.1853 


0.0160 


0.5027 


-0.1469 


0.1869 


0.1106 


-0.1188 




0.0713 


-0.1645 




-0.3165 


0.3084 




0.3207 


-0.1863 




-0.1699 


-0.1402 




0.0299 


0.0244 


0.3153 


0.1633 


-0.0703 


0.3874 


0.1037 


0.0021 


0.3874 


-0.1037 


-0.0021 


0.2121 


-0.4561 


0.3351 


0.0467 


0.0737 


0.0553 


0.4692 


0.2924 




0.3109 


0.4884 





experimental values is presented. 



B. CdTe in rock-salt structure 

In the second step we analyze the cadmium telluride in 
RS structure. To our knowledge, there are only few pa- 
pers, which report obtaining RS CdTe by applying high 
pressure to the ZB CdTe crystals. In Refs. |35| and [36| 
the optical studies of the band structure of this material 
are reported. The transition from ZB to the RS phase 



TABLE II. Experimental and theoretical values of the effec- 
tive masses of PbTe at the L point, mf, ml denote the lon- 
gitudinal and transversal effective masses of holes (h) and 
electrons (e), respectively. 





m* h /mo 


m$ h /mo 


ml J m 


m$ e /m 


Calc* 


0.162 


0.033 


0.133 


0.0281 


Calc3 


0.0799 


0.0133 


0.0799 


0.0107 


This work 


0.294 


0.0276 


0.272 


0.0241 


Expt.^2 


0.31±0.05 


0.022±0.03 


0.24±0.05 


0.024±0.03 




FIG. 4. The calculated within our tight-binding model band 
structure of PbTe along the symmetry lines of the Brillouin 
zone. 



is observed at (3.8 ±0.2) GPa as a dramatic decrease of 
the sample transmittance (the samples become virtually 
opaque between 3.9 GPa and 4.5 GPa). The measure- 
ments revealed that the change from ZB CdTe to RS 
CdTe results in a characteristic shift of the valence band 
maximum away from the T point towards the L and K 
points, due to p-d hybridization effects. As a result, two 
types of indirect gaps occur. However, it is very diffi- 
cult to determine experimentally the indirect band gaps 
of RS CdTe, because during the phase transition from 
tetrahedral to octahedral coordination the large number 
of defects and dislocations induced by the applied pres- 
sure form band-tails states.— 

The band structure of CdTe in RS structure was al- 
ready determined by using DFT-LDA method. In the 
calculation performed in Ref. |3fJ, RS-CdTe turned out 
to be a semimetal, in which the conduction band mini- 
mum at the X point would be nearly 2.5 cV below the 
valence band maximum located at S, midway the T - 
K line. This calculation, as well as another one pre- 
sented in Ref. [3|| does not take into account the spin- 
orbit effects, which are important in the analysis of the 
electronic properties of CdTe. The calculations of band 
structure with spin-orbit coupling were performed, us- 
ing linear muffin-tin orbital method, by Christensen and 
Christensen^l The band structure obtained in the lat- 



ter exhibits also a metallic character. Additionally, the 
band gap at the T point, at the center of Brillouin zone, 
is very close to zero. It should be emphasized that the 
authors of all the mentioned above results attribute the 
obtained semi-metallic character of RS-CdTe to the lim- 
itations of DFT-LDA methods only, i.e., to the fact that 
these calculations underestimate considerably the band 
gaps and can even give false overlaps between valence 
and conduction bands. Indeed, reflectance experiments 
in the mid-infrared do not indicate a metallic behavior of 
RS CdTe^ Giider et al. after a careful analysis of the 
experimental results obtained for RS-CdTe at different 
pressures made an assumption that RS-CdTe is a narrow 
gap semiconductor, with the energy gap of few hundreds 
of meV at the T point.— This assumption is consistent 
also with results of Ref. |35|. 

As described above, there are neither exact experimen- 
tal data nor a reliable theoretically calculated band struc- 
ture, to which we can fit our tight-binding parameters 
for RS CdTe bulk crystal. In this situation we decided 
to try to obtain a better description of the band struc- 
ture of RS CdTe in high symmetry directions of the fee 
Brillouin zone by using the ab initio procedure described 
in Sec. II. In this calculations the spin-orbit interactions 
are taken into account in the VASP code. The results are 
presented in Fig. [5j As one can notice in Fig. [5j we also 
obtain that the conduction band minimum of RS CdTe 
is located at the X point and is below the valence band 
maximum near the L point. The band structure is quali- 
tatively similar to the structure presented in Ref. |371 but 
in our case the band gap at T point equals ~ 0.805 eV. 
Thus, the presented in Fig. [5] structure while repeating 
the drawbacks of the other ab initio DFT structures, has 
the advantage of agreeing with the assumption of few 
hundreds of meV energy difference between the conduc- 
tion and valence bands at the T point. The parameters 
for this material presented in Tabic U are obtained by 
fitting the tight-binding band structure to the one pre- 
sented in Fig. 03 
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FIG. 5. RS-CdTe band structure along the main symmetry 
lines of the Brillouin zone, as calculated by using DFT-LDA 
method within the VASP code. 



C. MnTe 

The band structure of manganese salts, particu- 
larly MnTe, has been the subject of several different 
models.— ~— According to Allen et a/.,— the electronic 
structure of MnTe is determined by the combined effect of 
an exchange splitting of the Mn 3d states and of a strong 
hybridization of these states with the anionic p states. 
The hybridization tends to delocalize the d electrons and 
mediates the contribution of the d states to the structure 
of the valence and conduction bands. The impact of hy- 
bridization effects on the total band structure depends 
on the position of the d states. A full tight-binding de- 
scription of the electronic structure of RS MnTe, which 
is based on the spin-fluctuation theory, was presented 
by Masek et at— The obtained band structure exhibits 
metallic character with the Fermi level fixed within the 
half-filled d band. The top of the valence band is at 
L and the bottom of the conduction band is at X. In 
this tight-binding calculations, however, the spin-orbit 
interactions were not included and only the interactions 
between the nearest-neighbors were taken into account. 
This is not consistent with our description of PbTe bulk 
crystal, where we consider not only the interactions be- 
tween the nearest-neighbors but also next-nearest neigh- 
bor tight-binding integrals. For a proper description of 
Pbi-a-Mn^Te mixed crystals, in particular for applying 
the virtual-crystal approximation, it is most reasonable 
to take into account the same number of neighbors in 
both constituent materials. Thus, we decided not to 
take for the tight-binding description of MnTe crystal the 
paramctrization of Rcf. ,4lJ. Instead, we have considered 
a model with the nearest-neighbor cation-anion as well as 
next nearest anion-anion and cation-cation inter-atomic 
couplings, like done before for PbTe. Still, the model 
parameters were fitted to reproduce the band structure 
obtained by Masek et al^ The electronic structure of 
RS MnTe resulting from our model is shown in Fig. [6] It 
should be noted, that in our tight-binding description of 
RS MnTe the spin-orbit splittings could not be included, 
because they were not present in the structure in Ref. |4U 
which we tried to reproduce. The sp 3 d 5 tight-binding pa- 
rameters for RS MnTe are presented in Table |U 



D. Pbi-^CdzTe and Pbi-^Mn^Te mixed crystals 

To calculate the band structures of Pbi_ x Cd x Te and 
Pbi-a-Mn^Te alloys within the tight-binding approach, 
we apply the virtual crystal approximation, i.e., all the 
on-site energies and interaction integrals are assumed to 
be in x part equal to the given parameter for CdTe or 
MnTe, and in the remaining (1 — x) part equal to that 
of PbTe. In Fig. [7] the dependence of the energy gap 
of Pbi-^Cd^Te crystals in L, S and A points of Bril- 
louin zone on the Cd content, up to x = 0.2, is pre- 
sented. The calculated band gap in the L point is almost 
a linear function of composition and increases with x like 




FIG. 6. The calculated within tight-binding model (with 
nearest- and next-nearest-neighbors interactions taken into 
account) band structure of MnTe in RS structure. The red 
dashed line is the Fermi level. 



dEg/dx « 2 eV. It should be noted that this is ca 30% 
faster increase than that suggested by the experimental 
results of Ref. [lj- Fig. [8] presents the results of similar 
calculation for Pbi-^Mn^Te alloy for Mn concentrations 
also up to x — 0.2. As shown in the Figure, the energy 
gap in the L point increases with the Mn content ca 35 
meV/at %. The obtained increase of the band gap of 
Pbi-^Mn^Te with the Mn content is also slightly more 
rapid than the value 25 meV/(at % of Mn) estimated 
from the experimental data in Ref.[8|. On the other hand, 
however, the DFT calculations presented in Ref. [28| lead 
to probably too small slope of 15 meV/(at % of Mn). 
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7. (Color online) Dependence of the energy gaps in 
iCdxTe on the Cd content x. The blue squares, black 

stars and red dots denote the energy gaps at A, E and L, 

respectively. 



It should be recalled here that our study of PbTe-based 
mixed crystals was mainly motivated by the idea that in 
these materials a change of relative position of the heavy 
and light hole valence bands can lead to higher thermo- 
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FIG. 8. (Color online) Dependence of the energy gaps in 
Pbi-a;Mna:Te on the Mn content x. The blue squares, black 
stars and red dots denote the energy gaps at A, E and L, 
respectively. 
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electric power as compared to the p-type PbTe. It can 
lead to an increased ZT because the total thermoelectric 
power factor for the material derives from the contri- 
butions from all extrema . 12 ' 42 Also in Ref. M the strong 
increase of thermoelectric power found in Pbi-zMn^Te 
was explained by assuming a change in sign of the separa- 
tion energy between the band extremum of light holes at 
the L point and the band of heavy holes at S. A similar 
behavior can be expected in the Pbi_ x Cd x Te crystals. 
Indeed, in Figs [7] and [8] we observe that the energy dis- 
tance between the valence and conduction bands in the 
other extrema (at A and £) increases with adding Cd 
or Mn much slower than that in the L maximum. To 
study this behavior more carefully, we determine the en- 
ergy differences between the valence band maximum at 
the L point and another maximum at the S, as well as 
between L and A, i.e., we determine the so called "side 
energy gaps" . The obtained side gaps as a function of the 
composition x for Pbi-^Cd^Te are presented in Fig. [5] 
and for Pbi-^Mn^Te in Fig. [TO] 

First we notice that according to our calculations in 
Pbi_a;Mii x Te, although the side energy gap S - L dimin- 
ishes with x, up to x=0.2 the maximum at £ is lower 
in energy than the maximum at the L point of the Bril- 
louin zone. Thus, this is a slower decrease than that 
suggested in Refs |8| and [28| (within the ab initio calcu- 
lations with reduced spin-orbit interactions, presented in 
\2&l , it was obtained that for Mn concentrations higher 
than x ~ 0.09 the £ appears above the L maximum). For 
Pbi-^Cd^Te we obtain that for x > 0.1 the band maxi- 
mum at A appears above the top of the valence band at 
£ and that both side energy gaps should change sign for 
x = 0.2, as shown in Fig. |H1 Despite these quantitative 
differences, our results for both materials, Pbi-^Mn^Te 
and Pbi_ x Cd x Te, qualitatively agree with the idea that 
in these mixed crystals higher x enhances the role of the 
heavy holes from the vicinity of the other maxima of 
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valence band. This is especially valid for higher concen- 
trations of the p-type carriers in the samples. 



IV. SUMMARY AND DISCUSSION 

Our analysis of the stability of Pbi-^Cd^Te, per- 
formed by ab initio methods, shows that up to x ~ 0.8 
the total energy of the ternary alloy is lower for the RS 
structure. This result is consistent with the observation 
that it is easier to obtain rock-salt Pbi-^Cd^Te crys- 
tals with considerable amount of Cd than zinc-blende 
CdTe crystals doped with Pb. The lattice parameters 
of Pbi_ x Cd x Te and (Pbi_ x Mn x Te) mixed crystals de- 
crease with Cd (Mn) content x. The lattice parame- 
ters ao diminish with x< 0.12 like dag/dx = —0.434 A 



(Pbi.^Cd^Tc) and da /dx = -0.683 A (Pbi_ x Mn x Te), 
in good agreement with the experimental findings. 

The energy structures of Pbi-^Cd^Te and 
Pbi-^Mn^Te alloys have been analyzed using tight- 
binding description of constituent materials and virtual 
crystal approximation. For this purpose we have per- 
formed a tight-binding parametrization of rock-salt 
PbTc, CdTe and MnTc. In contrast to previous ap- 
proaches, our tight-binding parameters of lead tclluridc 
lead not only to correct overall band structure of 
PbTe in the whole Brillouin zone but also accurately 
reproduce the experimental bulk effective masses. Due 
to the lack of exact experimental data for RS CdTe, the 
tight-binding description of this material is based on the 
results of DFT calculations - still, our model recovers 
the predicted band gap in the center of the Brillouin 
zone. An increase of the L-point band gap with x has 
been obtained for both studied ternary alloys. The 
calculated band gaps in the L maximum are almost 
linear functions of the composition x and compare 
well with the experimental results. We have calculated 
also the energy differences between the valence band 
maximum at the L point and another maximum at S, 
as well as between L and A. In both studied materials 
the side energy gaps S - L and A - L diminish with x. 
Thus, for higher concentrations of the p-type carriers in 
the samples an enhanced contribution of the heavy holes 
from the valence band in the vicinity of the other valence 
band extrema is predicted. To discuss how adding Cd 
or Mn to PbTc changes the thermoelectric properties of 
the material, we have calculated the Seebeck coefficient 
(S) in both ternary alloys, using above described band 
structures of the mixed crystals. 

As a first approximation we consider Seebeck coeffi- 
cient within Mahan-Sofo theory, which is given by the 
simplified Mott expression^ 3 - 



7T 2 h 2 
3 Q 



1 dD(E) 
D dE 



1 dii(E) 1 



/' 



dE 



(1) 



E — Ef 



where Ef is the Fermi energy, q is the carrier charge; 
D(E) and fi(E) are the density of states and the mobility, 
both energy dependent. 

In Fig. [TT] we show how the first term in the Mott 
equation (which describes the role played by the changes 
of the density of states in the thermoelectric power) de- 
pends on the content of Cd in Pbi-^Cd^Te. The cal- 
culations have been performed for three different hole 
concentrations: p = 1 x 10 19 cm -3 , p = 3 x 10 19 cm -3 
and p = 1 X 10 20 cm~ 3 , i.e., different positions of the 
Fermi level in the valence band. As one can see in the 
Figure, the D~ 1 dD(E)/dE value increases dramatically 
whenever the top of the heavy hole band at the A point 
reaches the Fermi level. For p = 1 x 10 19 cm -3 this hap- 
pens for x ~ 0.16 and for p = 3 x 10 19 cm -3 for x w 0.12. 
For the highest, p = 1 x 10 20 cm -3 , hole concentration 
the heavy hole band with the top at £ contributes to 
D~ 1 dD(E)/dE value at the Fermi level even for x = 0. 



For this hole concentration the small increase at x = 0.06 
is attributed to the onset of the contribution from the A 
maximum (compare Figs [11] and [5]) . Here it is worth 
to recall that the concept of carrier pocket engineer- 
ing to produce convergence of symmetrically incquivalent 
bands has been suggested first for low-dimensional ther- 
moelectric nanostructures^ and then extended to bulk 
materials^ In these papers it has been suggested that 
the convergence of many charge carrying valleys has an 
effect of producing large m* without explicitly reducing 
the mobility fi and that a large valley degeneracy should 
improve thermoelectric properties of the materials. In 
Ref. |45| this effect was shown in PbTci-^Sea,, where the 
L and £ valence bands can be converged, giving an in- 
creased valley degeneracy of 16. We note that accord- 
ing to our calculations, in Pbi-^Cd^Te, under proper 
conditions, the valley degeneracy as high as 22 can be 
achieved, due to additional contribution from the sec- 
ondary valence band (A - with the degeneracy equal to 
6), which in this material is also very close to the L and 
£ bands. 
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The presented above Seebeck coefficient enhancement 
through density of states modification is a promising 
route ; 46 ' 47 but this approach risks the reduction of carrier 
mobility4£ The optimal electronic performance of a ther- 
moelectric semiconductor depends strongly also on the 
weighted mobility— To study Seebeck coefficient more 
precisely, in particular to take into account the mobil- 
ity term, we have followed the scheme proposed in the 
paper of Madsen and Singh^S i.e., we have performed 
the calculations using Boltzmann transport theory£2r— 
and constant relaxation time approximation^ 3 - The lat- 
ter approximation is based on an assumption that the 
scattering time determining the electrical conductivity 
does not vary strongly with energy. The advantage of 
using constant relaxation time is that the thermoelectric 
power can be directly calculated from the band structure 



as a function of carrier concentration and temperature, 
with no adjustable parameters. In Ref. |49| the thermo- 
electric transport tensors (i.e., the electrical conductivity 
<t(T, /x) and Seebeck coefficient S(T, /i)) are defined by 
the following expressions: 



Sa = (O 



ai^aj 



where 



<T a p{T,n) = ^ 



J af3 



(e)[- 



dU(T,e) 
de 



(2) 



}de (3) 



v a p(T,[t) 



eTVL 



<J afs (e)(e - n)[- 



i)s 



]de (4) 



and the transport distribution function tensor <7 a f3(e) is 
given by 

e 2 >-^ 

o a ti{e) = — 22Tv a {i,k)vp(j,'k)S(e - e iM ) (5) 



j,k 



where £i : k is band energy, v a (i, k) is the band velocity 
{dei^/dha); r, (j,, and /„ are the relaxation time, chem- 
ical potential, and Fermi-Dirac distribution function, re- 
spectively. More details about this model can be found 
in Ref. Ha 

The dependence of the thermoelectric power in p- 
type Pbi-^X^Te mixed crystals (X denotes either Cd 
or Mn) on the composition x was calculated for the 
same three carrier concentrations (p = 1 x 10 19 cm~ 3 , 
p = 3 x 10 19 cm -3 and p = lx 10 20 cm -3 ) and room tem- 
perature 300K. Eigenvalues of 200 x 200 x 200 k-points 
mesh were calculated by tight-binding method. The ob- 
tained results for Pbi-^Cd^Te are shown in Fig. [12] and 
for Pbi-^Mn^Te in Fig. Q21 As one can see in the Fig- 
ures, in both materials the Seebeck coefficient increases 
with increasing Cd (Mn) content. The thermoelectric 
power for p-type alloys with high carrier density differs 
appreciably from the corresponding values for PbTc. It 
should be noticed that this happens not only because of 
the greater significance of heavy-mass hole conduction, 
which results from the smaller energy separation between 
light and heavy hole valence bands. The increase results 
also from the more parabolic nature of the light-mass 
valence band (a consequence of the larger direct energy 
gap), as described already by Rogers and Crocker many 
years ago^ This is especially seen in Fig. Q2] because 
according to our models in Pbi-^Mn^Te the admixture 
of the heavy holes contribution is not as strong as ex- 
pected and it does not play a considerable role for the 
Mn content shown in the Figure. Still, according to our 
calculations, adding ca 10% of Mn to PbTc should en- 
hance the thermoelectric power by ca 15%. In PbTe with 
the same amount of Cd a much higher thermoelectric 
power should be observed. The results presented in Figs 
IT31 and [12] recover also the well known decrease of the 



Seebeck coefficient with the concentration of holes t 54 i 55 
which is due to the strong interdependence between dif- 
ferent properties of the thermoelectric material via the 
carrier concentration. 
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FIG. 13. (Color online) The dependence of the calculated 
thermoelectric power of p-type Pbi-^Mn^Te crystals on the 
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